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The Projected Configuration Interaction (PCI) method starts from a collection of mean-field wave 
functions, and builds up correlated wave functions of good symmetry. It relies on the Generator 
Coordinator Method (GCM) techniques, but it improves the past approaches by a very efficient 
method of selecting the basis states. We use the same realistic Hamiltonians and model spaces as 
the Configuration Interaction (CI) method, and compare the results with the full CI calculations in 
the sd and pf shell. Examples of 24 Mg, 28 Si, 48 Cr, 52 Fe and 56 Ni are discussed. 



O 
O 

(N 
C 

a 
in 



o 

=5 

> 
in 

CN 

o 

oo 
o 



X 



PACS numbers: 21.60.Cs,21.60.Ev,21.10.-k 



I. INTRODUCTION 

The full configuration interaction (CI) method using a 
spherical single particle (s.p.) basis and realistic Hamilto- 
nians has been very successful in describing various prop- 
erties of the low-lying states in light and medium nuclei. 
The realistic Hamiltonians, such as the USD [l|, Qin the 
sd shell, the KB3 @, FPD6 Q and GXPF1 [Jin the pf 
shell, have provided a very good base to study various 
nuclear structure problems microscopically. 

Some sd and pf nuclei such as 28 Si and 48 Cr are well 
deformed. Their collective behavior is confirmed by the 
strong collective E2 transitions and the rotational behav- 
ior of the yrast state energies, 75(7) ~ 7(7 + 1). The 
mean-field description in the intrinsic frame naturally 
takes advantage of the spontaneous symmetry breaking. 
This approach provides some physical insight, but the 
loss of good angular momentum of the mean-field wave 
functions makes the comparison with the experimental 
data difficult. The CI calculations in spherical basis pro- 
vide the description in the laboratory frame. The angular 
momentum is conserved, but the physical insight associ- 
ated with the existence of an intrinsic state is lost. 

There is a long-lasting effort to connect the mean-field 
and CI techniques. Elliott was the first to point out the 
advantage of a deformed intrinsic many body basis and 
developed the SU(3) Shell Model [| for sd nuclei. In El- 
liott's model, the rotational motion was associated with 
strict SU(3) symmetry, approximately realized only near 
20 Ne and 24 Mg. This model was limited only to the sd 
shell. 

The Projected Shell Model(PSM) Q can be consid- 
ered as a natural extension of the SU(3) shell model to 
heavier systems. In this model, the quadrupole force, 
that Elliott's model used, the monopole pairing and 
the quadrupole pairing forces were included in the PSM 
Hamiltonian. The deformed intrinsic Nilsson+BCS ba- 
sis are projected onto good angular momentum, and the 
PSM Hamiltonian is diagonalized in the space spanned 
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by the projected states. The Nilsson model Q has been 
proven to be very successful in describing the deformed 
intrinsic single particle states, and the quadrupole force 
was found to be essential for describing the rotational 
motion [f|. Despite its simplicity PSM was proven to 
be a very efficient method in analyzing the phenom- 
ena associated with the rotational states, especially the 
high spin states, not only for axial quadrupole deforma- 
tion, but also for the octupole Q and triaxial shapes 
[Tol . IllT l . However, its predictive power is limited because 
the mulitpole-multipole plus pairing Hamiltonian has to 
be tuned to a specific class of states, rather than an re- 
gion of the nuclear chart. 

Besides the Projected Shell Model, another sophisti- 
cated approach based on the projection method is MON- 
STER and the family of VAMPIRs [3. The MON- 
STER is similar to the PSM, but the basis is obtained 
by projecting the Hartree-Fock-Bogoliubov (HFB) vac- 
uum and the related 2-quasiparticle configurations onto 
good quantum numbers, including neutron and proton 
number, parity, and the angular momentum. VAMPIR, 
which performs the energy variation after the projection, 
is more sophisticated than MONSTER. However, the 
particle number plus angular momentum projection im- 
plemented in MONSTER and VAMPIR requires at least 
3-dimensional integration in the axial symmetric case. 
This makes it very difficult to extend these models to 
non-axial cases, where the projection with 5-dimensional 
integration would be needed. MONSTER and VAMPIR 
can use realistic Hamiltonians. 

Instead of using quasiparticle configurations, as PSM 
and MONSTER/ VAMPIR do, the Quantum Monte 
Carlo Diagonalization (QMCD) method [13| takes the ad- 
vantage of the Hartree-Fock (HF) mean-field that breaks 
the symmetries of Hamiltonian. The mean-field is not 
restricted to axial symmetry. QMCD starts from an ap- 
propriate initial state, as in the shell-model Monte Carlo 
approach [lj], to select an optimal set of basis states, 
and the full Hamiltonian is diagonalized in this basis. 
More basis states are iteratively added until convergence 
is achieved. It is very interesting that only the M- 
projection is applied to in the basis states, yet the total 
angular momentum seems to be fully recovered from the 
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diagonalization when the convergence limit is achieved. 
However, the process of selecting the basis states is re- 
portedly very time consuming. 

The Generator Coordinate Method (GCM)[lj| is a 
standard method to describe collective states that goes 
beyond the mean-field. The angular momentum projec- 
tion technique is a special case of GCM when the genera- 
tor coordinates include the Euler angles. In the standard 
GCM, the excited states are constructed in a basis of 
HFB vacua, which differ in one or several collective coor- 
dinates. However, in many cases, it has been found that 
the excitation energies provide by GCM are too high. 
The GCM wave functions may be more appropriate to 
describe more correlations in the ground state, rather 
than the excited states. However, our investigations (see 
below) show that neither the g.s. nor the excited states 
are accurately described by this simple GCM procedure 
when realistic effective Hamiltonians and HF vacua are 
used. 

In the present work, we developed a new method called 
the Projected Configuration Interaction (PCI), which 
uses GCM and PSM techniques. The PCI basis includes 
angular momentum projected states generated from a 
class of constraint HF vacua. However, in addition to 
these GCM-like states PCI includes a large number of 
low-lying np — nh excitations built on those states. In 
PCI the deformed intrinsic states are Slater Determinants 
(SD), and the particle number projection is no longer re- 
quired. This method can use the same realistic Hamil- 
tonian as any CI method. This feature makes the direct 
comparison between CI and PCI possible. One advan- 
tage of this method is that it keeps the number of basis 
states small, even for cases when the CI dimensions are 
too large for the today's computers. 

The paper is organized as follow. Section II presents 
the formalism used by PCI. The efficiency of the angular 
momentum projection is discussed in section III. Section 
IV analyzes the choice of the PCI basis and presents some 
results for the sd and pf nuclei. Section V is devoted to 
the analysis of the quadrupole moments and E2 transi- 
tions. Conclusions and outlook are given in section VI. 



II. THE METHOD OF THE PROJECTED 
CONFIGURATION INTERACTION (PCI) 

The model Hamiltonian used in CI calculations can be 
written as: 



(1) 



i>j,k>l 



where, cj and Cj are creation and annihilation operators 
of the spherical harmonic oscillator, et and Vijki are one- 
body and two-body matrix elements that can be obtained 
from effective interaction theory, such as G-Matrix plus 
core polarization [16j ; it can be further refined using the 
experimental data[5|, [TtJ ■ 



One can introduce the deformed single particle (s.p.) 
basis, which can be obtained from a constraint HF solu- 
tion. Alternatively, one can take the single particle states 
obtained from the following s.p. Hamiltonian, 
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Here 



h sp h — ^ E, 



c]ci 



(2) 



(3) 



is the spherical single particle part of the Hamiltonian 
with the same eigenfunctions as the spherical harmonic 
oscillator, but the energy of each orbit is chosen such 
that the result from this Hamiltonian is near to the HF 
solution; and £4 are the quadrupole and hexadecupole 
deformation parameters, p = -J 
and we take [H, fl9( | 



m^-r is dimensionless, 



hiu = 45A 1/3 - 25A 2/3 



(4) 



The deformed s.p. creation operator is given by the 
following transformation: 



(5) 



where the matrix elements Wki = (bk\ci) are real in 
our calculation. Inserting the reversed transformation 
of Eq. © 



4 = E^ & £ 



(6) 



into the model Hamiltonian Eq. |T]), we obtain the trans- 
formed H in the deformed s.p. basis 



h?l,b]b%b k . 



X] U ijkl u i"j 
i>j,k>l 



(7) 



Here and are one-body and two-body matrix 
elements of H in the deformed basis. The Slater Deter- 
minant (SD) built on the deformed single particle states 
is given by 



\K) = \s,e) = blbt 2 . 



(8) 



where s refers to the Nilsson configuration, indicating the 
pattern of the occupied orbits, and e is the deformation 
determined by 62 and €4. The matrix element (k\H\k'), 
can be calculated using Eq.([7]). 

Generally \k) doesn't have good angular momentum 
/. It is well known that the projection on good angular 
momentum would add correlations to the wave function. 
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The general form of the nuclear wave function is therefore 
a linear combination of the projected SDs (PSDs), 



\^im) — ^ IikkPmk i K ) ; 



(9) 



Kk 



where 



f«K = ^ / >nw' Uh ai)im) no) 

is the angular momentum projection operator. D\ [K \s 
the D-function, defined as 



D'mkW = (IM\R(Sl)\IK)*, 



(11) 



R is the rotation operator, and is the solid angle. If one 
keeps the axial symmetry in the deformed basis, D I MK in 
Eq lgnU reduces to d{ IK ((3) = (IM\e~^ J y\IK) and the 
three dimensional fl reduces to (3. 

The energies and the wave functions [given in terms of 
the coefficients Jj Kk in Eq.©] are obtained by solving 
the following eigenvalue equation: 



^ {Hk K ,K>k' - E I N K K: K'K')flKK' - 0, (12) 
K'k' 

where Hr k k'k' ano - ^Kk K'k' are ^ ne matrix elements of 
the Hamiltonian and of the norm, respectively 

Hkk,k'k' = WHPkx'W), (13) 
NL,K',' = Mk'W)- (14) 

The matrix element of Hj iK K , K , and K , K , can be 
expanded as 



2/+ 1 



dOI}j fK ,(n)<re|-H'.R(Q)|re'), (15) 



8vr 2 
( k \Pkk'\ k ') 

f d^ KK ,{si)(K\h{n)\K'), (is) 



where 

(k\HR(Q)\k') = ^2(k\H\k")(k"\R(Q)\k'). (17) 

k" 

Here k" has the same deformation as re, and runs over all 
SDs with nonzero (k\H\k"). Therefore, the main problem 
is to calculate the matrix element of the rotation R(Cl) 
between different SDs, {k\R(Q)\k') . 

| re) and \k') may have different shapes. We denote 
a\a the single particle operators with deformation e a 
that create the SD |re), and w ,b, with deformation q,, 
that create |re'). Expressing {k\R(CI)\k') with a\a and 



b\ b, one gets 

(re|-R(0)|re') 
= (|o iB ...o ia o il A(n)6t i 6t j ... 6 tj) 
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Here 



= £w^(Q)4 

i V Ik ) 



(19) 



Note that the rotation never mixes the neutron and pro- 
ton wave functions. Therefore, the matrix elements in 
Eq.(fr5]) can be calculated for neutron and proton sepa- 
rately. 

The choice of the PCI basis will be described in Section 
IV. 



III. THE EFFICIENCY OF THE PROJECTION 
ON GOOD ANGULAR MOMENTUM 

Let's analyze the quality of the deformed single particle 
states for the USD Hamiltonian Q ■ The £7, energies (see 
Eq.Q) are properly adjusted so that the SD with the 
lowest energy is close to the HF vacuum, as will be dis- 
cussed in this Section. For sd shell, we use Ed 5/2 = —3.9 



MeV, E s = -1.2 MeV and E d 



1.6 MeV, and the 



corresponding Nilsson levels as functions of quadrupole 
deformation £2(64 = 0) are plotted in FigfTJ 

To find the lowest value i?def of (k\H\k), the energy 
surface of 



E exp (s,e) = (s,e\H\s,e) 



(20) 



needs to be calculated for each configuration s. Here, we 
use the USD Hamiltonian, and take 28 Si as a case study. 
The SD with all 12 particles occupying the Nilsson orbits 
originating from the d 5 / 2 spherical 
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FIG. 1: The Nilsson levels for sd shell 



-125 



-130 




FIG. 2: (Color online) Energy surfaces of the USD Hamilto- 
nian with respect to the unprojected SD(so) (upper surface) 
and the projected SD(so) with 1=0 (lower surface), as func- 
tions of quadrupole £2 and hexadecupole €4 deformation. 



orbit are used in Eq. (20). Such an SD is denoted 
by so hereafter. This configuration provides the lowest 
energy for a large range of deformations and is expected 
to describe very well the mean-field minimum, i.e. the 
ground state (g.s.), (k\H\k) at certain deformation(s). 
The energy surface of the so configuration as a function 
of (£2,64) is shown as the upper surface in Fig. [2j The 
minimum energy, E^d — —129.55 MeV at deformation 
(e2,£4)dcf = (—0.44,-0.10), which is very close to the 
HF energy (E HF = -129.61 MeV) [Hj]. Besides, the cal- 
culated quadrupole moment, — 70.1fm 2 , is also very close 
the the HF value, — 69.6fm 2 [l8|. Thus, we reached an 
approximated solution to the HF mean-field from a sim- 
ple Nilsson Hamiltonian. For other deformed sd nuclei 
the situation is similar to 28 Si, as shown in Table HI 

The ground state of 28 Si should have good angular mo- 
mentum with spin 1 = 0. We calculate the expectation 
value of H with respect to the projected SD: 



(s,e\HP^ K \s,e) 



AP*kk\ 



(21) 



where K is defined by J z |s, e) = K\s, e). The energy sur- 
face of E oxp (I = 0, Srj,e) is shown as the lower surface 
in Fig. H The minimum value E pi (I = 0) = -133.64 
MeV at (e 2 ,e4) P j = (-0.40,-0.20), which is about 4 
MeV lower than the HF energy, and 2.30 MeV above 
the exact ground state energy of the USD Hamiltonian 
at —135.94 MeV. Similar calculations for other deformed 
sd-shell nuclei are reported in Table. |TJ These results 
indicate that the restoration of the angular momentum 
has a significant contribution to the ground state (g.s.) 
correlation energy. 

An even stronger argument in favor of imposing good 
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FIG. 3: The ground state energies (upper panel) and the 
corresponding (J 2 ) values (lower panel) of 28 Si from spherical 
and deformed CI as functions of different np-nh truncation 
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TABLE I: Energies (in MeV) with USD Hamiltonian for some deformed sd Nuclei. -E sp h is the spherical HF energy, -Ehp is the 
minimum of the deformed HF energy, Edei is the lowest energy of Eq. (|20|l and (e%, e.4)def is the deformation for -Edef, E p j(I = 0) 
is the lowest energy of Eq. (|21fl at spin I — and (e2, £4) P j is the deformation for E p j(I — 0). -Efuii ci is the exact solution of 
the USD Hamiltonian. 



Nucleus 


T7i a 
-fisph 


771 a 
£/HF 




(C2, £4)dcf 


£ pj (I = 0) 


( e 2, £4)pj 


-EWn ci 


2U Ne 


-31.79" 


-36.38 


-36.38 


(0.60, -0.04) 


-39.86 


(0.56, -0.26) 


-40.49 


24 Mg 


-68.20 


-80.17 


-79.80 


(0.50,0.08) 


-83.04 


(0.58,0.04) 


-87.09 


28 Si 


-126.03 


-129.61 


-129.55 


(-0.44,-0.10) 


-133.64 


(-0.40, -0.20) 


-135.94 


36 Ar 


-222.75 


-226.56 


-226.56 


(-0.32,0.08) 


-229.78 


(-0.20,0.30) 


-230.51 



"Data taken from Ref-Ejj. 

"The -25.46 MeV in Ref. p]J has been corrected as -31.79 MeV. 



angular momentum to the wave functions can be found 
by comparing the results of CI calculations using spher- 
ical and deformed s.p. bases, as shown in Fig. [3] The 
deformation for the deformed CI basis is (e2,£4)dcf = 
(-0.44,-0.10) as shown in Table Q] The energies at Op- 
Oh are the values of E sp h and E^ef shown in the same 
Table Q] The E de { is much lower than E sph , but (J 2 ) 
for the deformed SD is as large as 20.35, far away from 
(J 2 )=0 of the spherical HF state. At the 2p-2h trun- 
cation level, the spherical CI energy drops very close to 
that of the deformed CI, while the latter is struggling 
for recovering the angular momentum, and (J 2 ) drops 
to 11.30. At the 4p-4h truncation level, the spherical 
CI energy is even lower than the deformed ones, yet the 
angular momentum of the latter hasn't been completely 
recovered. The advantage of the deformed mean field is 
then completely lost in the CI due to the lack of good 
angular momentum. At the 8p-8h truncation level the 
deformed CI state completely recovers its angular mo- 
mentum, and the results from spherical CI and deformed 
CI become equivalent. Thus, we conclude that the de- 
formed CI won't benefit from the deformed mean field 
unless the angular momentum is recovered from the out- 
set. The PCI method described in Section II maintains 
the advantages of the deformed mean field and of config- 
uration mixing specific to CI, and could provide further 
improvement to the E p] (I = 0) in Table |TJ 

To describe those non-zero spin states, one can change 
the spin from to I in Eq. (21), and keep s and e un- 
changed. Then one gets a sequence of energies, E p j(I) 
shown as D=l in Fig. [U One can see that the values of 
Epj(I) have approximately reproduced the /(/ + 1) fea- 
ture of a rotational band in the low spin region. The 2 
MeV gap between the D=l band and the Full CI band 
is expected to be reduced by taking advantage of the 
CI techniques, in which the interaction between differ- 
ent intrinsic configurations has been properly considered. 
Results of the PCI for different bases are also shown 
in Fig. |U First, we take 41 projected SDs having the 
same configuration (sq), but different shapes, as in the 
GCM method. The 41 shapes considered are taken from 
€2 = —0.6 to 0.6 by a step of 0.03, and £4 is obtained 
by minimizing the energy of Eq. (|21[) with / = at each 
(-2- One can see that the calculated g.s. band is not 
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FIG. 4: Calculations of the ground state band in 28 Si with 
different basis. Except for the GCM:D=41 case, the deforma- 
tion for all the other cases is (ei = —0.40, £4 = —0.20) 



improved. This can be understood by the large overlap 
between the PSDs in the GCM basis, therefore, the effec- 
tive dimension of the basis is far below 41. Adding more 
PSDs with new shapes in the basis doesn't make much 
improvement. Therefore, we considered up to 2p — 2h 
excitations of the deformed sq configuration used in the 
D=l case, and built a PCI basis by selecting D=40 SDs 
(see next section). With this basis we solved the gen- 
eralized eigenvalue problem described in Section II. The 
results (PH:D=40) in Fig. 0] indicate a significant im- 
provement over GCM (GCM:D=41). As the number of 
excited PSDs increases up to 107, we obtain the band 
(labeled by PH:D=108), which is only about 300 keV 
above the band obtained with full CI. Therefore, we con- 
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FIG. 5: Comparison of the calculated ground state bands in 
28 Si with deformations (£2,£4) p j = (—0.40,-0.20) [used in 
FigB] and (e a , £4)def = (-0.44, -0.10), as shown in TableU 



elude that the PCI method is not only an efficient way 
of reducing the dimension of the CI calculations, at least 
in deformed nuclei, but can be also very instrumental in 
making the connection between intrinsic states and lab- 
oratory wave functions. 



IV. CHOICE OF THE PCI BASIS 

As already mentioned, one important problem is the 
choice of the PCI basis. It is very difficult to find a 
set of SDs that would make an efficient PCI basis, due 
to the complexity of possible structures that exist in the 
nuclei, such as the spherical states, rotational states, var- 
ious vibrational states, shape coexistence, etc. It would 
be helpful to be able to dentify the most relevant intrinsic 
states of a nucleus when we start a PCI calculation. Un- 
fortunately, there is a large number of ways of achieving 
that goal. Below, we describe some of them. Assuming 
that we found these N starting states that we denote as 
\kj, 0) (j = 1, ...A), we further consider for each selected 
\hij, 0) a set of relative np-nh SDs \hij, i), and select some 
of them into the basis according to how much effect they 
can have on the energy of the state. Therefore, the gen- 
eral structure of the PCI basis is 



' Op — Oh, np — nh 
|ki,0), 

|*t2,0), \K2,ii),--- 

|kjv,0), \n N ,ii),-- 



(22) 




-0.6 -0.3 0.0 0.3 0.6 

e 2 



-0.6 -0.3 0.0 0.3 0.6 

e 2 



FIG. 6: (Color online) £ cxp (s,e) values in 24 Mg and 28 Si, as 
functions of quadrupole deformation £2. Each curve repre- 
sents a Nilsson Configuration s, while £4 was chosen to min- 
imize _B exp (s,e) for each £2- The absolute minimum for each 
s is indicated by a black dot. 



A is somewhat arbitrary, but it can be increased until 
reasonable convergence is achieved. Some examples are 
given below. As a preliminary application of PCI, we con- 
sider the same set of SDs \kj, 0) in Ed.(|22p for all angular 
momenta /. The deformation of each configuration s may 
be chosen from the minima of either Eq. (I20|l or Eq. (j2l"j) 
at certain /. However, there seems to be not much differ- 
ence between these two choices, especially if the studied 
nucleus is well deformed. As shown in Fig. \5\ the PCI 
results of 28 Si with deformations (ea, £4) P j and (e2, C4)def, 
(see Table IJ) are quite close to each other. Therefore, for 
well deformed nuclei it would be simpler to use Eq. (j2"0|) 
for determining the shape of each starting configuration 
Sj. In the present work we only consider axially sym- 
metric shapes. Some low-lying curves of E axp (s, e) in Eq. 
(|2U)) for 24 Mg and 28 Si are plotted in Fig. [5] as functions 
of 62. £4 was chosen by minimizing the E exp (s, e) for each 
£2- Each configuration Sj corresponds to a curve in this 
figure. The minima are marked by dots. 

In order to describe the low-lying states, we first choose 
those SDs whose energies are minima of E Qxp (s, e), shown 
as dots in Fig. [SJ and take them as starting 0) states 
in Eq. (|22p. Secondly, we try to include CI- like correla- 
tions. Each 0) may have correlations with its lp-lh 
and 2p-2h SDs, through the non-diagonal matrix 

elements of H. For each j we include those SDs that 
satisfy the following criterion: 

1 



AE = -(E 



Ei + v / (^o-^) 2 +4|F| 2 ) > E cuU (23) 



where E Q — (k,0\H\k,Q), Ei — (n,i\H\n,i) and V = 
(K,0\H\n,i). (We skipped the subscript j to keep nota- 
tions short.) 
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FIG. 7: Energies for 24 Mg and 28 Si calculated using PCI (open 
circles) and full CI(filled circles). 20 energies for each spin are 
presented. 



In this way we constructed the basis for the PCI calcu- 
lations of 24 Mg and 28 Si in the sd shell. For the 24 Mg, we 
included TV = 40 lowest energy \nj,Q) SDs with \K\ < 2, 
and 1749 SDs are selected for the basis when E cut = 10 
keV. For the 28 Si case, N = 60, \K\ < 4 and E cut = 10 
keV selects 3366 SDs. The choice of \K\ < 2 seems to be 
arbitrary, and at this point it is related to the number 
of states one wants to include in the basis. One should 
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FIG. 8: The lowest 5 energies at each spin for 48 Cr calculated 
using PCI (open circles) and full CI (filled circles). 
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FIG. 9: Calculated yrast band energies using PCI (open cir- 
cles) and full CI (filled circles) for 52 Fe and 56 Ni. The de- 
formed band in 56 Ni starting around 5 MeV is also shown 
(see text for details). 



emphasize that the PCI method combines the advantage 
of using a set of s.p. bases of different deformations with 
that of the CI np—nh configuration mixing. In m-schcmc 
CI, for example, one can get the states for all spins us- 
ing the lowest M value (0 for even-even cases), although 
using higher M values could be somewhat more efficient 
for higher spins. In PCI one has to weigh in the advan- 
tages of using a larger number of np-nh excitations with 
those of using larger K values for larger spins. Using our 
choices for the K values the lowest 20 PCI energies of 
each spin are compared in Fig. [7]with the full CI results. 
Almost all low-lying energies for each spin reproduce the 
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full CI results very well. We conclude that the selected 
PCI basis has already carried almost the whole informa- 
tion of the low-lying states, yet we chose only a small 
number of SDs as compared with those used in the full 
m-scheme CI. 

The QMCD [l3| method also uses a relatively small 
number of selected SDs. For a basis of 800 SDs QMCD 
found a g.s. energy of —89.91 MeV for 24 Mg. In the 
present calculation, after the orthogonalization of the 524 
K = SDs, 453 SDs are used to calculate the 1=0 en- 
ergies (see Table [TTT]) . The resulting PCI g.s. energy is 
-89.94 MeV. Furthermore, the QMCD paper only re- 
ports the 6 lowest states with spin up to I = 4, while 
our PCI calculation with the selected 1749 SDs can pro- 
vide tens of excited states with spin running from to 12 
that compare well with the full CI results. Our approach 
seems to largely extended the range of nuclear states that 
can be calculated using a relatively small basis. 

The results for 28 Si are similar to those of 24 Mg. Com- 
paring with the N = 1 calculation, that has been dis- 
cussed in Section III, the dimension is much larger be- 
cause one needs to describe a large number of low-lying 
states of the same spin. For 1 = the number of K = 
\nj, 0} SDs is N = 15, and a total of 747 independant 
SDs are selected from the 886 K = SDs. The first 
20 + states, up to about 18 MeV excitation energy, are 
reasonably close to the full CI results. For I ^ the 
situation is very much the same. 

The pf calculations show similar trends. For 48 Cr, the 
CI m-scheme dimension is 1963461. Full CI calculations 
can reproduce the rotational nature of the yrast band 
and its backbending [2(J. With PCI we choose N = 30, 
\K\ < 4, and for E cut = 0.2 keV 5973 SDs were selected. 
The lowest 5 energies for each spin are compared with the 
full CI results in Fig. [5] Again, we reached a satisfactory 
agreement between PCI and full CI. 

52 Fe and 56 Ni are not good rotors, and the present 
choice of the basis is not so efficient as for 48 Cr. How- 
ever, one can still select a PCI basis that can provide a 
reasonable description of the yrast energies when com- 
pared with the full CI calculations, as shown in FigJ5] 
The number of deformed SDs used in these calculations 
is about 7000, which is just a small fraction of the full 
CI m-scheme dimension for 52 Fe, and for 56 Ni (see Table 
III). However, the choice of the PCI basis for 52 Fe and 
56 Ni needs to be improved. Details will be presented in 
the forthcoming paper. 

For 56 Ni a rotational band has been identified in exper- 
iment [2l| that starts at around 5 MeV excitation energy. 
The band was successfully described by full CI calcula- 
tions in the p/-shell using the GXPF1A interaction (T^. 
It was also shown [17j that this band has a 4p-4h char- 
acter, but in order to get a good description of the band 
lOp-lOh excitations from (/7/2) 16 configuration are neces- 
sary. In the present PCI calculation this rotational band 
has been identified (see Fig. 9) in the calculated excited 
states by inspecting the quadrupole moments and the 
B(E2) transitions (See Table HI)) . The results seem to 



TABLE II: Quadrupole moments (in efm 2 ) and B(E2,I — > 
I — 2) (in e 2 fm 4 ) values of the states in the deformed band of 
56 Ni calculated with PCI and CI. The energies for 1^0 are 
relative to the 7 = bandhead. The CI values are taken from 
Ref. H3| 



Spin 


Energy (MeV) 


Q(i) 


B{E2) 




PCI 


CI 


PCI 


CI 


PCI CI 





-200.295 


-200.762 










2 


0.425 


0.395 


-48.9 - 


41.6 


567.1 413.2 


4 


1.312 


1.241 


-44.0 - 


55.2 


646.0 598.0 


6 


2.541 


2.448 


-70.5 - 


56.2 


708.0 609.3 


8 


4.103 


3.991 


-74.5 - 


47.2 


938.6 558.4 



TABLE III: PCI Dimensions for several spins compared with 
those of full CI for the nuclei described in this Section. 



Nucleus 




1 = 


1 = 2 


1 = 4 


1 = 6 


24 Mg 


SM 


1161 


4518 


4734 


2799 




PCI 


453 


1569 


1579 


1492 


28 Si 


SM 


3372 


13562 


15089 


9900 




PCI 


747 


2266 


2919 


2711 


48 Cr 


SM 


41355 


182421 


246979 


226259 




PCI 


1671 


3888 


5618 


5594 


52 Fe 


SM 


1.8 x 10 6 


8.0 x 10 6 


1.2 x 10 7 


1.2 x 10 7 




PCI 


4288 


6422 


6811 


6806 


56 Ni 


SM 


1.5 x 10 7 


7.1 x 10 7 


1.1 x 10 8 


1.1 x 10 s 




PCI 


2373 


3968 


5452 


6231 



be in reasonable good agreement with those of the full 
CI calculations. The details about the calculation of the 
quadrupole moments and the B(E2) transitions are given 
in the next section. 

The number of SDs selected for the PCI calculation 
are shown in the corresponding figures and compared 
with the full CI m-scheme dimensions. However, the ef- 
ficiency of the PCI truncation can also be assessed by 
comparing the PCI dimensions with the coupled-I CI di- 
mensions. In PCI the states used for the diagonalization 
of each spin in Eq. (12) are obtained from the selected 
PSDs after the latter are orthogonalized and the redun- 
dant states are filtered out. The remaining states form 
the PCI basis for each spin. These dimensions are com- 
pared in Table ITTT1 with the full coupled-I CI dimensions. 
One can see that PCI can provide reasonable good results 
using much smaller dimensions than the corresponding 
m-scheme and/or coupled-I CI calculations. 

Although the PCI calculations could have much 
smaller dimensions than those of the corresponding full 
CI calculations, the H and N matrix in Eq. (12) are 
dense. To get an idea of the workload of a PCI cal- 
culation, the time necessary to calculate the H and N 
matrices for 56 Ni at 1=0, 2, 4, 6, 8 takes around 8 hours 
using one processor, and the calculation of 5 eigenvalues 
for each spin takes about 2 hour, when the generalized 
Lanczos method is used. As a comparison, the full CI cal- 
culations of the same states using the modern coupled-I 
code NuShellX 22] could take several weeks, when one 
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processor is used. One should also observe that the PCI 
computational load is shifted towards the calculation of 
the matrix elements, which can be more efficiently par- 
allelized than any large matrix eigenvalue solver. 



V. QUADRUPOLE MOMENTS AND BE2 
TRANSITIONS 



The quadrupole moments and B(E2) transitions are 
jiven by 



I6tt 



Q{I) = \/-g-(*/M=l|Q20|*JAf=j) 



I6tt 



— (II, 20|J7)<* / ||Q 2 ||* J ) (24) 
5 



30 
20 

a ■ 

S-io 

o 

-20 
-30 



24 Mg(USD) 



I 1 I ' I ■ I 



50 
40 
•I 30 h 
20 
10 



28 Si(USD) n 48 Cr(GXPF1A) 



• Full CI 
o PCI 
• . • 




2 4 6 8 1012 2 4 6 8 1012 4 8 12 16 
200 200 400 




2 4 6 8 1012 

Spin 



2 4 6 8 1012 
Spin 



4 8 12 16 

Spin 



B(E2;I i ^I { ) = -L^-\(^ If \\Q a \\^ 

= IM&pz,)! 2 (25) 

where 

Qa m = e„ (p) r 2 F v (0,0), (26) 

with effective charges taken as e„ = 0.5e, e p = 1.5e. 
The reduced matrix element in Eq. ([24"f and (|25| can be 
expressed in terms of the PCI wave functions 

<*/ / IIQ 2 ||*/ ! ) 

= ^ fliKuflfK'n'RlfK'K'JiKu, (27) 
Kk,K'k' 

where 

RlfK'n'JiKu 

V 

Explicit expression for the ($ K / \Qx u Pj^t_ v k \$k) can be 
found in Ref. 0]. Here (I t K' - v,Xv\I f K') is a Clebsh- 
Gordon coefficient. 

Using the PCI wave functions for the states shown in 
FigEland Fig0 we calculated the Q(I) and the B(E2) 
values for the yrast band in 24 Mg, 28 Si and 48 Cr. The 
results are shown in Fig llOl The very good agreement 
between PCI and full CI indicates that the PCI wave 
functions are almost equivalent to the exact ones. 

VI. CONCLUSIONS AND OUTLOOK 



FIG. 10: Calculated quadrupole moments (upper panel) and 
B(E2) (lower panel) from PCI (open circles) and full CI (filled 
circles)for 24 Mg, 28 Si and 48 Cr. 

that of the spherical basis, due to the loss of rotational 
invariance that is not recovering fast enough. 

Therefore, we propose a new strategy called Projected 
Configuration Interaction (PCI) that marries features of 
GCM, by projecting on good angular momentum SDs 
built with deformed s.p. orbitals, with CI techniques, 
by mixing a large number of appropriate deformed con- 
figurations and their np — nh excitations. This ansatz 
seems to be very successful in describing not only the 
yrast band, but also a large class of low-lying non-yrast 
states. This method seems to work extremely well for 
deformed nuclei, and it needs further improvements to 
the choice of the PCI basis for spherical nuclei. 

For the sd and p/-shell nuclei that we studied we show 
that the quadrupole moments and the BE(2) transition 
probabilities are very well reproduced, indicating that 
not only the energies, but also the wave functions are ac- 
curately described. This suggests that this method could 
be helpful in finding relations between the intrinsic states, 
which offer more physics insight, and the laboratory wave 
functions provided by CI. 

One of the simplifications that makes these calcula- 
tions with tens of thousands of projected states possible 
is the approximation of the deformed mean-field with an 
axially symmetric rotor defined by its quadrupole and 
hexadecupole deformations. This approach seems to de- 
scribe very well the HF mean field for the sd and p/-shell 
nuclei. We plan to extend our method to include octupole 
deformation, and to consider triaxial shapes. 



In this paper we investigate the adequacy of using a de- 
formed s.p. basis for CI calculations. We show that even 
if the starting energy of the deformed mean-field is lower 
than the spherical one, the series of np — nh truncations 
in the deformed basis exhibits slower convergence than 
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